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Abstract 

We study the time required for a wave packet to tunnel beyond a square 
barrier, or to be reflected, by envisaging a physical clock which ticks only 
when the particle is within the barrier region. The clock consists in a magnetic 
moment initially aligned with the x axis which in the barrier region precesses 
around a constant magnetic field aligned with the z axis, the motion being in 
the y direction. The values of the x and y components of the magnetic moment 
beyond or in front of the barrier allow to assign a tunneling or reflection time 
to every fraction of the packet which emerges from the barrier and to calculate 
tunneling times tt, x and 7~T,y and reflection times tr jX and TR iy . The times tt, x 
and r^y (tr x and TR jy ) are remarkably equal, and independent of the initial 
position (in front of the barrier) of the packet. 



PACS numbers: 03.65.Bz,73.40.Gk 



I. INTRODUCTION 



An old problem in Quantum Mechanics which has attracted the attention of several 
authors (see M- j2|- |3j for extensive reviews) is that of a satisfactory definition of the 
time required for such a peculiarly quantum mechanical phenomenon as tunnelling. In a 
previous paper || we investigated the behaviour of a wave packet impinging onto a square 
barrier and proposed definitions of the interaction time (or dwell time) td, as well as of the 
transmission time tt and the reflection time tr which were based on the time dependence 
of the probabilities Pi (i = 1,2,3) of finding the particles in the region x < —d, \x\ < d, 
x > d respectively, the potential barrier lying between x = —d and x = d. 
According to these definitions 
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each of the times td, tt and tr are weighted sums of the time intervals dt. For td the 
weight is the probability of finding the particle within the barrier, for tt and tr the weight 
is the fraction of the total transmitted (reflected) packet which at time t has not yet been 
transmitted (reflected), T and R being the probability of the particle being definitely 
transmissetted and reflected at t — > oo respectively td is independent of the choice of the 
lower extremum t 2 in eq. |T], whereas tt and tr depend on the values of t\ and £3, which 
should be chosen as the time when the packet begins to interact with the barrier. In ref. 
[|J we fixed t\ and £3 as the time r t such that the contribution to the dwell time of the 
interval [0,r e ] is equal to that fraction e of the total dwell time which can be regarded as 
the resolution of our procedure. We estimated this resolution as O.OItd. 

It is desirable, however, to remove also the slight arbitrariness implicit in the above 
choice, and to have a more objective definition of the times involved in the tunneling 
process. To this purpose, we resort to the well known device (see refs. 0- 0- 0) of 
attaching to the particle a clock which runs only when the particle is in the barrier 
region. This clock can be provided by a magnetic moment which precesses around a 
constant magnetic field B confined within the barrier region. In order for this interaction 
not to affect substantially the barrier height, it is necessary to consider the limit of 
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infinitesimal magnetic field, which will cause an infinitesimal precession. The measure 
of this precession, i.e. of the average of the spin of the particle in the region beyond the 
barrier (or the average over the reflected component) should yield the time the particle 
has spent in the barrier before being transmitted (or reflected). So far, the use of a 
physical clock for determining the tunneling time has been limited to the case of stationary 
solutions. 

In order to work out this program, we have studied the evolution of a wave packet 
representing a particle with spin 1/2, polarized in the x direction, initially located in front 
of the barrier, which travels in the y direction towards the barrier. The average kinetic 
energy is lower than the barrier height. In the barrier region the particle interacts with 
a magnetic field B directed in the z direction. As a consequence, the transmitted and 
reflected packets will be no longer polarized in the x direction. 

A magnetic moment fx initially aligned in the x direction which emerges at time t 
from the barrier region after interacting a time r(t) with a constant magnetic field has a y 
component \i y = — /xsin(u;i) where (in units with h — 1) to — fiB is the Larmor frequency. 
In turn, the x component is fi x (t) = /j,cos(out). We extract r(t) from the average, at time 
t, of the operators /i x and /i y over the transmitted packet, which are proportional to the 
average values over the transmitted packet of S x and S y , which we denote as (S x (t))3 and 
(S y (t))s. We solve for r(t) either of the equations 
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i / sin[cjr(a;)]— —dx 

The above equations are to be viewed as the definition of r(t). The tunneling time is 
calculated as follows: 



/ r(t) 
Jo 



dP > it 



In other words, we weight every increment dP 3 of the probability of having the particle 
beyond the barrier with the time r(t) that same fraction dP% of the transmitted packet has 
spent within the barrier, the time being measured through the spin precession. Through- 
out the calculation we take advantage of the infinitesimal value of u, so that we are content 
with the first non- vanishing term. 
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We notice that in principle eqs. |] and [| yield different definitions of r(t), i.e. we have 
two different observables, S x and S y , whose measurement allows an indirect measurement 
of r(t). Thus, eq. ^ could yield two different results for the tunneling time tt. As a 
matter of fact, the results obtained starting from eq. |] or eq. |5] are equal well within the 
accuracy of the calculation, although, as will be discussed in sec. 4, the values of r(t) are 
not exactly equal (but they are practically equal starting from sufficiently early times). 
Moreover, the results are independent from the location of the packet at time t = 0. 

With analogous considerations we derive the reflection time tr, starting from the 
average of S x and S y in the region in front of the barrier: 

r * , dP x 
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where the label 1 denotes the average over the region in front of the barrier. Also in this 
case the values obtained by choosing either of the values of r(t) extracted from eqs. 
and |H 

f°° , .dP 1 , 

_ _ Jo at ( Q \ 

" P 1 (co) (9) 

are very near and do not depend on the initial position of the wave packet. 

In conclusion, the prescriptions embodied in eqs ^-^| (0-§3|) appear to give a rather 
objective measurement of the tunneling time tt (of the reflection time tr). The resulting 
values are in agreement with those we found in ref. Q and do not suffer of the ambiguity 
in the choice of the lower extremum in eqs. |2| and [| i.e. the time when the particle begins 
to interact with the potential. 

In sec. 2 we present the problem, and in sec. 3 we derive the time r(t) and the 
consequent values for tr and rx- We comment our results in sec. 4, while sec. 5 is 
devoted to the conclusions. 



II. DEFINITION OF THE PROBLEM 

We consider a wave packet describing a particle of mass m, moving in the y direction, 
where a potential barrier lies between y = —d and y = d: 
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V(y) = V 



k 2 
2m 



\y\ < d 



(10) 



A homogeneous magnetic field B directed in the z direction is confined in the barrier 
region. The particle has a magnetic moment fi = gS which in the barrier region is coupled 
to the magnetic field by the interaction fi ■ B = gBS z . Consequently, the hamiltonian of 
the particle is (we use units with % = 1): 
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The initial state of the particle is described by the wave function 



4>(y,o) = <f>(y)x 

where x is the spinor (l/y / 2)(l, 1), such that S x \ '- 
located in front of the barrier and peaked at y = y Q . 
function whose evolution was studied in ref. HI, 



(12) 



<P{y) = / a(k)ip k (y)dk 



- x an d (j>(y) is a Gaussian packet 
We choose for 4>(y) the same wave 

(13) 



where ipk(y) are the eigenfunctions of the hamiltonian |TT] with u = and 

25^ 1/4 
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The parameters involved in the problem are chosen as follows: 

1, d = V2, K v = 9.9, d=2, k = 10, y = -15 



m 



(14) 
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In order to study the time evolution of the wave function [L2| under the action of the 
hamiltonian [11] we must expand 12 in the basis of the eigenfunctions ipk,+ {y) and ipk,-(y) 
with S z = 1/2 and S z = —1/2 respectively. These functions are the eigenfunctions 
of the square barrier hamiltonian, where the barrier height is Vq — uj/2 and Vq + uj/2 
respectively. We note however that the initial wave function is located in the region 
y < —d, where the dependence of ipk,+(y) and ipk-iy) on S z is only through the reflection 
coefficient A(k). The contribution of this part of the wave function to the coefficients 
a + {k) = J ipl + (y)cj)(y)dy and a_(/c) = / tpl _(y)(p(y)dy is negligible, due to the Gaussian 
shape of 0. Consequently, we can safely assume the coefficients a + (k) and a-(k) to be 



the same, equal to a(k) as given in eq. [14. 
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The time evolution of the wave function [12| is given by 



_ a{k)i, K+ {y)e^ k2t ' 2m dk 
that is 

^(y, t) = ^J a(k)My)e- lk2t/2m dk (16) 



III. EVALUATION OF THE TUNNELING TIME AND THE REFLECTION 

TIME 

We derive the tunneling and reflection times according to the program outlined in the 
Introduction. We focus our attention on the calculation of the tunneling time, the case 
of the reflection time being quite similar. 

Let -Ps(t) be the probability of finding the particle beyond the barrier (y > d) at time 
t (0 is the Heavyside function): 

P 3 (t) = (i;(y,t),e(y-d)i;(y,t)) (17) 

We want to associate to the increment dP^, = dP^/dtdt of this probability the time r(t) 
which this fraction of the transmitted packet has spent in the barrier region. The tunneling 
time for the packet will be calculated according to the equation || 
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The time r(t) will be extracted from the average of the spin component S x or S y at time 
t in the region beyond the barrier: 

, s = (il>(y,t),Q(y - d)S x ip(y,t)) ^ 



(jj(y,t),e (y-d) S y ^(y,t)) 
Ps(t) 



(S y (t)) 3 = WWh^-yvVW)) (19) 



The key for deriving r(t) is the spin precession which occurs when the particle is in the 
barrier region where a constant magnetic field aligned with the z axis is located. For a 
magnetic moment fi aligned with the x axis at time t — and interacting with the field, 
the time dependence of its x- and y-component is given by 
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H x — ficos(ut) and fi y = — fism(u)t) (20) 

where u = [iB. The same time dependence holds for the x and y components of the spin 
angular momentum: 

1 1 

S x = -cos(ut) and S y = — - sm(uit) (21) 

If the fraction dP 3 of the packet emerging beyond the barrier in the time interval dt at 
time t has spent a time r(t) in the barrier region, we assign to this fraction a contribution 
to (S x (t)) 3 and (S y (t)) 3 equal to 

1 dP 

dS x = -— -cos[ujr(t)}dt (22) 
1 dP 

dSy = ~2~lt sin [ cjr ( t )] rft ( 23 ) 
Consequently, the values of (S x (t)) 3 and (S y (t)) 3 will be given as 

^ = 2^15m (24) 

< 5 »«> 3 = -2^W < 25 » 



We have to extract r(t) from either of eqs. [24] and |25|. This can be done by multiplying 



both sides of the equations by ^(t) and taking the derivative with respect to t. It is useful 
to observe, however, that we need to consider the limit when uj tends to zero, since we 
are interested in the tunneling time beyond a barrier of heigth Vo, and a non-infinitesimal 
magnetic field would alter the barrier heigth. 
If we write 

{Sx{t)h= W) and {Sy{t)h= W) (26) 



we get: 



dN x (t) 1 dP 3 (t) 

^r = 2— cos[ujT{t)] (27) 

dN y (t) _ 1 dP 3 (t) 



dt 2 dt 



sm[ur(t)} (28) 
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In eq. ^ we approximate cos(u)t) ~ 1 — \u) 2 t 2 , and in eq. ^ we put sin(u;t) ~ cut. We 
note also that N x and P3 are even in u, whereas N y is odd in u. Moreover, at order in 
uj we have N x = ~P 3 , N y = 0. By putting 

N x = iV (0) + oo 2 NP 



P 3 = Ni 0) + cu 2 P. 



2p(2) 



from eqs. ^ and ^B] we get 
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dt 



dt 2 



r(t) 2 



dt 



(29) 



dNW(t) ldN^(t) 



r(t) 



(30) 



dt 2 

According to the above equations, we have two different possible prescriptions to derive 
the time r(t) which must be inserted into eq. We denote by tt, x the result obtained 
from eq. |6] using eq. |29|, which yields for r(t) 



T(t) 
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and by r Ty the result for r(t) obtained from eq. |30|: 

T(t) = 



dt 
dNW 



dt 

The calculation of tt )X requires the calculation of the integral 
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whereas, according to eq. || and eq. [30|, the value of TT, y depends only on the values at 
t = 00 of iV"W and P 3 , and is given immediately by 
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Performing the calculations, we find: 

tt x = 3.53 



dt 



r T,- 



P,foc) 



3.52 



(34) 



(35) 



In order to find tr we repeat the above calculations starting from the values of 



(S x (t)h 



N x (t) 
Pi(t) 



and 



(Sy(t)h 



Ny(t) 
Pl(t) 



where now 



N x (t) = (iP(y,t),Q(-y-d)S x i>(y,t)) 



and 



NJt) = Mv,t)M-V-d)S y il>(y,t)) 



These allow the extraction of r(t) according to formulae analogous to eqs. ^T] and p2| , 
where of course P3 is substituted with Pi. The reflection times tr x and tr iX are given by 
the equations 
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r R,y = 

The results are as follows: 



Pi (00) Jo 



dt 



^(00) 
Pi (OO) 



TR, X = 0.52 



r Riy = 0.56 



(37) 



(3f 



Finally, we observe that, while eqs. |35] and |38] directly connect the average value at t = 00 
of S y in the region in front or beyond the barrier with tr iX or tt, x; this connection does not 
hold for the x-component, i.e. we do not have the relation, which holds in the stationary 
problem M, 



1 



- 2 2 



(39) 



As a consequence, there is no use in introducing t z as t z = 2(S z )/u, since the condition 
Tg = rl + Ty, which is derived from the identity (S x ) 2 + (S y ) 2 + (S z ) 2 = 1/4, no longer 



holds, due to the failure of eq. 39. 
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IV. COMMENTS 



The results obtained in the previous section for both the tunneling and the reflection 
time using r x (t) or r y (t) are near to each other, and are remarkably similar to the values 
found in ref. ||, where we found rx = 3.39, tr = 0.55. This supports the reliability of 
the definition of the tunneling and reflection times proposed therein and the conclusions 
about the impossibility of extending the notion of tunneling (reflection) time proposed 
for stationary problems by several authors to the case of wave packets. For a comparative 
discussion of some results of this kind we defer the reader to ref. [|J . 

Some comments are in order with regard to the dependence of the tunneling and 
reflection time on the initial position of the packet and on the different ways to evaluate 
this time (i.e. the choice between r Tx and r Ty , or between r Rx and r Ry ). 

Concerning the first point, it is easy to see that the integral in eq. ^ can be written 

as 
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where D(k) is the transmission coefficient and the derivative with respect to uj is taken 
at uj = 0. No dependence on the value yo where the packet is peaked at t = is left. 

The comparison between tt jX and 7r, y (t RjX and r Riy ) is more delicate. We will focus 
our discussion on the tunneling time, the considerations for the reflection time being 
quite similar. Eqs. ^7] and |28| show that the condition for the identity of the values of 
r(t) derived from them is 
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This condition is not verified for any value of t. However, the ratio 
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(42) 



is less than 0.1 starting from t ~ 2 (less than 0.01 starting from t ~ 2.7), where the 
main contribution to the integral |G] comes from. This can be easily understood. With 
the parameters of the problem, the packet travels with a velocity approximately equal to 
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k& v /iri ~ 10. Given that the initial peak is located at y = —15, up to t ~ 1.5 the fraction 
of the packet in the barrier region is practically zero. 

As to the reasons why the LHS of eq. ^l] is very nearly equal to the RHS, while 
not being exactly equal, they can be traced to the following observations. Each term in 
the numerator of eq. ^ can be written in terms of the corresponding flux calculated at 
y = d. If ip + and ip' + ijp'_) are the values in y = d of the component with S z = 1/2 
(S z = —1/2) of the wave function and its derivative with respect to y, the numerator A" 
in eq. ^2] reads 

2 



N 



ip'*ip--r+ip-~-\j + +j-\ 2 (43) 



j + and j_ being the currents due to ip + and calculated in y — d. If the function 
ip + ijp~) is expressed in polar form, = p+e l9+ , (■?/>_ = p_e* 9_ ), we have j + = p 2 + Q' + 
(j- = p 2 _9'_) and A^ can be written as follows: 

N = - p +P '_f - (p^; 2 - p 2 J 2 ){ P \ - pi) (44) 

At order in u N vanishes, whereas the second order contribution in uo is due only to the 
first order variation, i.e. it can be obtained by writing p + = p — udp, p_ = p + udp and 
the likes for the other terms. As a consequence eq. can be written as follows: 



R = p ^ («) 

On the other hand, using eq. |27| and ^ and denoting by t x the value of r(t) in eq. 



and by r y the value of t in eq. EHl at second order in to R can be written as 



R = u\t 2 v - t 2 ) = 2u 2 r av Ar (46) 



y 

At being r y — t x and r av the average of t x and r y . As a consequence 
At 2 (p'dp - pdp') 2 - 8p 2 6'dp{6'dp + pdff) 



(47) 



'av P V 'av 

Since dp/ p and the analogous terms for p' appearing in eq. E7I are much smaller than 1, 
we have Ar/r av 1, and consequently tt, x = 7"T,y • 

V. CONCLUSIONS 

We have studied the process of a wave packet tunneling beyond a potential barrier 
which has a clock running only in the barrier region. This is achieved by the well known 
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device, so far envisaged only for stationary problems, of considering the particle endowed 
with a magnetic moment coupled to a magnetic field confined in the barrier region. The 
time is measured by means of the precession of the magnetic moment around the field 
direction. 

The particle is initially polarised in a direction orthogonal to the direction of motion 
and that of the field. So, for any fraction of the packet emerging beyond the barrier 
(or reflected by it), the value of each of the two components of the magnetic moment 
beyond the barrier (or in front of it) orthogonal to the field can yield a measurement of 
the time spent in the barrier region. We find that the results obtained using any of these 
components are quite similar, within the accuracy of the calculation. The results for the 
tunneling and reflection times obtained by means of this method are in good agreement 
with the results obtained for the same packet, based on the consideration of the time 
behaviour of the probabilities of having the particle in front or within or beyond the 
barrier Q. 

The tunneling and reflection times obtained in this way are independent of the position 
at time t = of the packet, provided it is located sufficiently far from the barrier edge. 
Thus, the definition of the tunneling or reflection time by means of a physical clock looks 
a promising approach when one wants to go beyond the monochromatic approximation. 
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